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A simple approach allowing to construct closed-form an- 
alytical zero-field magnetization distributions in cylindrical 
particles of a small thickness and an arbitrary shape (not 
necessarily circular) is presented. The approach is based on 
reduction of the non-linear Euler equations for magnetization 
vector field to the classical linear Riemann-Hilbert problem. 
The result contains all the distributions minimizing the ex- 
change energy functional and the surface magnetostatic con- 
tribution exactly, except for the neighbourhood of topological 
singularities on the cylinder faces where the result is approx- 
imate. The completeness of the analysis permitted to find 
a new type of a topological soliton in the case of circular 
cylinder. Also, an example of magnetic vortex in a triangular 
cylinder is given to investigate the role of the particle corners. 
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Small magnetic particles shaped as flat arbitrary cylin- 
ders with the size of the order of 100 nanometers and their 
arrays recently captured a lot of attention due to their 
unusual magnetic properties (as unique systems where 
two-dimensional topological solitons can be directly ob- 
served) and also because of their applications in the mag- 
netic random access memory devices. 

In large ferromagnetic bodies the equilibrium magne- 
tization distributions are usually represented as a set 
of magnetic domains separated by possibly bent one- 
dimensional domain walls 0|. Unfortunately, this ap- 
proach is hard to apply to small magnetic particles if 
their size is such that the domain wall does not fit in- 
side, resulting in considerable distortions of its profile 
with respect to the one-dimensional case. Consequently, 
it is necessary to consider this distribution as a whole. A 
number of finite-element simulations of the problem were 
already performed, see e.g. |^ and references therein. 

In this Letter an analytical approach allowing to 
find closed-form expressions for equilibrium topologically 
charged magnetization distributions in soft magnetic par- 
ticles shaped as flat arbitrary cylinders is presented. 

Mathematical formulation of the problem. Consider 
a particle made of an isotropic magnetic material, the 
anisotropy of soft material is neglected here. In the case 
of no applied magnetic field its total energy has two 
contributions, the exchange and the dipolar (magneto- 
static). Let us limit ourselves to the cylindrical particles 
with thickness L < Le = y^C/Mg (it usually can be 
relaxed |^ for L up to several Le), where C is the ex- 
change constant and Ms is the saturation magnetization 
of the material. In this case the magnetization distri- 



bution TO(r) — {mxTiTiY ^mz} ~ M{r)/Ms, |^i('^)| = 1 
can be assumed uniform along the cylinder thickness {Z 
axis), with r — {X,Y} the radius vector in the cylinder 
plane. In this notation the total energy of the particle (a 
functional of to) can be written in the following way 
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where V is the area of the cylinder face, V — 
{d / dX,d / dY,d / dZ} is the gradient operator, /iD[m] = 
H£i[fh]/Ms is the demagnetizing field (a functional of to) 
created by the magnetization distribution. The demag- 
netizing field can be expressed using the Maxwell equa- 
tions through its scalar potential Ho = WU{f), which 
is, in turn, a solution of the Poisson equation V^C/ = p 
with the requirement (due to the finite size of the parti- 
cle) that both \r\U and |r|^|V[/| are finite as |r| — > oo, 
p — W ■ M is the density of magnetic charges (on the 
boundary equal to the normal component of M) , Q . 

The exact solution for equilibrium TO(r) corresponds 
to the minimum of the functional (|^) and is given by 
a system of non-linear integral partial differential equa- 
tions. They are hard to handle even with the help of 
computer for a particular case. Thus, a simplification of 
the original problem (|l]) needs to be applied. 

To make such simplification let us qualitatively ana- 
lyze contributions of the different energy terms in (^ . It 
turns out that for sufficiently small particle there is a well 
defined hierarchy of energies. The most important term 
is the exchange, as when the particle gets smaller not 
only the amount of possible magnetic charges decreases 
(on par with the exchange energy) but also (due to ab- 
sence of magnetic monopoles) the positive self energy of 
charges gets more and more compensated by their neg- 
ative interaction energy. It is also clear, as the role of 
surface increases for small particle sizes, that the sur- 
face magnetostatic contribution is more important then 
the volume one. The surface contribution can be sub- 
divided into the energy of the magnetic charges on the 
faces and on the side of the cylinder. Because for flat 
cylinders the faces have usually larger area than the side, 
the former energy is more important. Thus, in flat cylin- 
ders we have the following energy hierarchy: exchange, 
face charges, side charges, volume charges. Note that the 
exact criteria for the necessary smallness of the cylinder 
strongly depend on its shape (due to the long-range char- 
acter of the dipolar interactions) and should be analyzed 
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on a case-by-case basis. Let us just assume for now that 
the cyhnder satisfies the criteria. 

The solution for the problem exactly minimizing the 
exchange energy functional ([^) with hjj = 0, and also 
having negligible amount of the face magnetic charges 
present only at the points of topological singularities was 
given in j|] as 
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where f{z) is an arbitrary analytic (in the sense of 
Cauchy-Riemann conditions [^) function of a complex 
variable, line over a symbol denotes complex conjugation, 
straight brackets denote absolute value of the complex 
number, z = X + iY is the complex coordinate, z = 1. 
The components of the magnetization vector are given by 
mx+imY — 2w/{l + \w\^) and ~ (1 — |wp)/(l-f |u)p). 
This solution consists of two parts both locally extremal 
to the exchange energy functional: the one convention- 
ally termed as soliton Q at \f{z)\ < 1 and the meron 
11. The role of the soliton "hat" (|/(z)| < 1) is to cover 
the topological singularity of the meron this is the 
only region of (^ where \mz\ > 0. Pure solitons are 
not realized in considered particles due to excessive face 
magnetic charges 

Let us now see if using the freedom of choice for f{z) 
in (^ we can exactly remove the side magnetic charges 
(whose energy is next in the hierarchy). This problem 
can be formulated as: to find a function f{z) analytical 
in the region V in such a way that Re [/(C)'^(C)] = (no 
magnetization components normal to the side), where C G 
C = dT> is the boundary ofD, and n{C) = n^iC) + iny{C) 
is the complex normal to C. Which makes it reduced to 
the well known linear Riemann-Hilbert problem . 

To write the solution let us denote the conformal trans- 
formation of the interior of the unit circle \t\ < Ito z € V 
as z = T(t), according to Riemann theorem such trans- 
formation exists for any simply connected region V. Un- 
der such transformation C transforms onto the boundary 
of the unit circle |A| = 1 and the normal is given by 
nxiX) = n(T(A)). Then, all the solutions for f{t) are 
given by 
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where oq, ai, 02 are arbitrary complex constants, 
F+{t) = F(t),\t\ < 1 and F-{t) = F{t),\t\ > 1 are 
the values of the Cauchy-type integral (^) having discon- 
tinuity at \t\ = 1. The presence of exactly three arbitrary 
constants is due to the index Q ind[— TiT/nT] = —2. 

The derivative of the conformal transformation T'{t) = 
G'(t)e*'^^*^ gives the scaling (the real function G{t)) and 



the rotation angle (the real function (p{t)) performed by 
a transformation at a point t, and A (when |A| = 1), is 
a normal to the unit circle. This allows to express (up 
to the real multiplier) the normal to the face boundary 
as nriX) ~ AT' (A). The real multipher in uxiX) is ir- 
relevant because only the combination nf/nT enters (Q). 
Then, the inte gral (|) gives F+(f) = -nr - logT'(i) and 
F-{t) = -logT'(l/t), where the identity 1/A = A and 
the residue theorem were used. Substititing these into 
(^) gives the final answer 



fit) = {itc + A- At^)T'{t), 



(5) 



where c = Im ai is a real constant and A = oq — 02 is 
a complex constant, the role of these constants is best 
illustrated by the following example. 

Circular cylinder. In this case T{t) = t and the mag- 
netization distribution is (0) with 
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The values of c and A need to be determined from the 
detailed analysis of the energy including magnetostatic 
and can not be determined based on the pole avoidance 
principle alone. However, the interesting point is that 
apart of the well known vortex solutions ||^ correspond- 
ing to |c| > 2\A\ there are also solutions of another type 
with |c| < 2\A\ corresponding to a pair of skyrmions |l^] 
(also known as hedgehogs) bound to the cylinder sides, 
see Fig. |l|. These solutions are relevant to the consider- 
ation of quasi-uniform magnetization states (by allowing 
the pair of skyrmions to move away from the cylinder 
boundary) and of a vortex nucleation. The nucleation of 
a vortex can be described as symmetrical moving of the 
pair of skyrmions from infinity (uniform state) to the par- 
ticle boundary with quasi-uniform states forming in the 
process and, then, by changing cx) > |A| > c/2 alongside 
the particle boundary to the point where two skyrmions 
annihilate and the vortex forms, moving subsequently to 
the particle center as A changes down to c/2 > \A\ > 0. 
This process will be considered in detail in a forthcom- 
ing paper. Let us now see what happens if a particle has 
corners by considering the next example. 

Triangular cylinder In this case the conformal map 
is given by the Schwartz-Christoffel integral T{t) — 
C /q(1 - t^y^/^dt. The side length of the triangle is 
1 if C = 22/3V3^/r(l/6)/r(l/3), r{x) is Euler's gamma 
function. Using (w we get 
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where T~^{z) means inverse function of T{t) which is 
unique because the transformation is conformal. This 
expression diverges at the corners of the triangle be- 
cause corners represent topological singularities of the 
type cxd/cxd in the meron. To avoid them we shall cover 
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corners by a soliton "hats" similarly as it was done H] 
for the singularities of the type 0/0 corresponding to the 
vortex centers. Thus, instead of (0) we shall use 
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J{Z)/C2 \f{z)\>C2, 
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where C2 > 1 is an additional arbitrary real constant. 
The result is plotted in Fig. |[ 

Summary. A general method is provided for obtain- 
ing closed-form analytical expressions for magnetization 
distributions in flat ferromagnetic cylinders of arbitrary 
shapes made of isotropic material. This problem is shown 
to be equivalent to finding the conformal transformation 
of the unit circle to the shape of the cylinder face T(t). 
Then, the magnetization distributions are given by (|5|) 
and (H) (or (^) if the particle has corners). No other mag- 
netization distributions exactly minimizing the exchange 
energy of the particle and having no surface magnetic 
charges on its side exist. A new type of soliton solu- 
tion was found by applying this procedure to the circular 
cylinder, see Fig. |l| (and also experimental Fig. 1(c) in 
[ prf , however interpreted differently there). 

The only term not included into this consideration is 
the magnetostatic energy of the volume charges but, as 
it was qualitatively argued, the effect of this term is the 
smallest for small particles and is expected to introduce 
only minor corrections to the magnetization distributions 
obtained here. If desired, these corrections can be found 
by perturbative expansions or by numerical analysis. 
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FIG. 1. The magnetization distribution (in-plane compo- 
nents) within the circular cylinder having two skyrmions 
bound on the side (||), c = 0, A = e'''^^/2. The dotted arrows 
show the solution continued into the region with no magnetic 
material to confirm that the magnetization non-uniformities 
are indeed skyrmions. 
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Magnetization distributions within triangular 
: a) with two side-bound skyrmions c = 0, 
= 4; b) with a vortex c = 1, ^ = 0, C2 = 4. 
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